Mon. Not. R. Astron. Soc. 000, 000-000 (0000) 



Printed 27 January 2012 



(MN IMftX style file v2.2) 



Infra-red emission line tomography of the intergalactic 
medium during the Epoch of Reionization 



(N 

o 

(N 
(N 



o 

u 

Oh; 

6 

o3 



Koki Kakiichi, Avery Meiksin, Eric Tittley 

SUPA* , Institute for Astronomy, University of Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, 



UK 



27 January 2012 



ABSTRACT 

The first major star-forming galaxies and Active Galactic Nuclei will produce Balmer 
and higher order extended haloes during the Epoch of Reionization through the scat- 
tering of Lyman resonance line photons off the surrounding neutral intergalactic gas. 
The optical depth dependence of the scattering rates will produce a signal sensitive to 
both the density and velocity fluctuations of the gas, offering the possibility of prob- 
ing the ionization region and flow field surrounding young star-forming galaxies. The 
requirements for detecting the haloes in the infra-red using a space-based telescope 
are discussed, along with an assessment of the possibility of detecting the haloes using 
the Tunable Filter Imager on the James Webb Space Telescope. 

Key words: atomic processess - cosmology: theory - line: formation - radiative 
transfer - infrared: general - scattering 
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1 INTRODUCTION 

The breaking of the redshift 2 = 7 barrier in the campaign 
to discover high redshift galaxies is closing in on the Epoch 
of Reionization (EoR), when sufficient stars, with a possi- 
ble contribution from Active Galactic Nuclei (AGN), formed 
to reionize the hydrogen produced in the Big Bang follow- 
ing the Recombination Era. In addition to the handful of 
sources with spectroscopically confirmed redshifts z > 7 
l|Vanzella et al.ll201ll : lOno et al.ll20ll ISchenker et~aT]|2012l ) , 
several plausible candidates have been identifi ed with higher 
photo metric redshifts extending up t o ■ z ^ 9 (McLur e et al.l 
l201ll h and possibly as high as z ~ 10 (|Bouwens et al.ll201ll V 
Measurements of the Cosmic Microwave Background (CMB) 
by the Wilkinson Microwave Anisotropy Probe ( WMAP) 
confine the EoR, if a sudden event, to z r — 10.4 ± 1.2 (lcr) 
|Larson et al.|[2oTlh . Thus the highest redshift galaxies dis- 
covered may have already entered the EoR. 

In addition to ionizing the IGM, high redshift galaxies 
will illuminate still neutral intergalactic hydrogen before the 
EoR has ended. The Lya radiation escaping high redshift 
sources will scatter in the vicinity of the sources, producing 
Lya hal oes with a characteris tic core radius of 70 kpc at 
z = 10 (|Loeb fc Rvbickilll999f) . These haloes would be ob- 
servable in the near infra-red, with characteristic diameters 
of 30 arcsecs. The sources will be overwhelmed by zodiacal 
light, but a dedicated space-based facility with narrow band 
imaging could plausibly detect a halo surrounding a very 
bright source. 
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Higher energy photons emitted by the continuum of the 
sources will scatter as well as they redshift into local Lyman 
resonance frequencies. Radiative cascades following the scat- 
tering of high order Lyman series photons will produce radi- 
ation from higher o rder hydrogen sequences like the Balmer 
and Paschen series {Meiksin 2 01oT ). The haloes will be highly 
extended because of the distances over which higher energy 
photons may travel before scattering. The Ha signal from 
within the inner 1 arcminute of the halo could, as for the 
Lya signal, be detected by an infra-red detector in space. 
The Balmer haloes, however, have the advantages over the 
Lya haloes of being both more extended and comprised of 
multiple orders. Combining the images for the different re- 
solvable orders could substantially reduce the required inte- 
gration times for a given source. 

Since the identical structures would give r ise to 

a 21cm signature arou nd a source |Madau et al.l 1 19971 : 
iPritchard fc Loebl l201lh , correlating the images with ra- 
dio detections using an instrument like the LOw Frequency 
ARray (LOFAR0 or the Square-Kilometre Array (SKa£| 
would further enhance the detections and probe the struc- 
ture of the underlying IGM as well as the ioni zation geom- 
etry produced by the source (|Tozzi et al.ll2000l h 

The purpose of this paper is to compute the expected 
Balmer signals allowing for the structure of the IGM. Since 
the scattering rate depends on the optical depth from the 
source to Lyman resonance line photons, the signal will 
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depend on the density, temperature and peculiar velocity 
structure of the IGM. The influence of inhomogeneities in 
these quantities is computed using a cosmological numerical 
simulation. 

The paper is structured as follows: the next section 
summarises the mechanism of the effect and the simulation 
model. A presentation of the resulting images follows. The 
observational requirements for detecting the signal are then 
discussed, followed by our conclusions. 



2 THE PRODUCTION OF BALMER HALOES 
2.1 The Lyn scattering rate 

The local Lyn radiation field within a gas parcel arises from 
two contributions, the direct UV photons emitted from the 
central radiating source redshifted to the local Lyn reso- 
nance frequency by the cosmological expansion and any ra- 
dial peculiar velocity v r of the gas, and the photons pro- 
duced in radiative cascades following the scattering of higher 
order Lyman series photons. The effect of redshifting limits 
the distance freely streaming photons from a source may 
travel before scattering. From the perspective of a gas par- 
cel at redshift z, only sources within the Lyman horizon 



i . horiz 
1 + Z n 



(1) 



are able to excite a hydrogen atom into an electronic stat e 
with principal quantum number n jBarkana fe Loebll2005h . 
Higher energy photons will have passed through a higher 
order Lyman resonance. As a result, the maximim possible 
distance of sources able to produce a direct scatter of Lyn 
photons decreases with increasing n. 

Higher order Lyman photons will scatter within the 
Doppler cor<jf]- The mean free path of Lyn photon within 
the Doppler core is 



n^ 2 A 



Bh(z)(7„ 

0.101(1 + z)' 3 T 1/2 n(n 2 - 1) pc 



(2) 



where hh(z) is the hydrogen number density, 
f n b/c is the Doppler width with Doppler parameter 
b = (2/cT/niH) 1 ^ 2 for gas temperat ure T, and a „ = 
(ire 2 /m e c)f in for oscillator strength /i n (|Meiksinl2 009). The 
approximation fi n ~ 1.56n~ 3 was used, accurate to bet- 
ter than 10 per cent, for n > 4. The mean free path is 
much smaller than the Jeans length Xj = ^(n/GpM) 1 ^ 2 — 
16(l + 2:) _3/ ' 2 r 1 ' /2 kpc for a sound speed c 3 and cosmic mass 
density pu over which the physical properties of the IGM 
will be nearly uniform. The redistribution of photon energy 
from scatters may thus be considered to occur locally, con- 
fined within a homogeneous medium. A higher order photon 
will scatter only a few times before decaying into a lower or- 
der photon, with a survival probability of ~ 0.8 per scatter 
l|Hiratall200rj ; iPritchard fc Furlanettol l2006L so that the ef- 
fects of spatial and frequency diffusion, redshifting and evo- 



3 At z = 8, for an IGM temperature T > 10 K (100 K), Ly<5 
(Ly7) and higher order Lyman resonance line photons scatter in 



the Doppler core iMeiksinll2uTol ). 



lution of the IGM will be negligible before the photon is 
destroyed. 

Including photons produced in cascades from higher or- 
ders, the net scattering rate of Lyn photons is 



Pn = 



pw + v v . p , 

n' >ra 



(3) 



for a maximum order n max of directly scattered Lyman res- 
onance line photons, where P n is the direct scattering rate 

P n 1} = Pt c (0)S n , (4) 

with P n nc (0) — a n L„ n (4-Kr 2 j hu n ) the scattering rate at the 
luminosity distance vl from a source of specific luminos- 
ity j^ii-fcagj/ fi+z) assuming no intergalactic attenuation 
(Meiksin 2010). Here, S„ accounts for the intergalactic at- 
tenuation and p n i >n is the probability that a Lyn' photon 
converts to a Lyn photon per scatter. 

The scattering rate may be expressed more succinctly 
in terms of the direct scattering rates as 



P n = 



(1) 



(5) 



where the scattering cascade matrix C n ' jn has been defined 
according to 



n 

c n ., n = C n ,, n „r, n „, B , 

n" >n 



(6) 



= for n > n! and r\ n ^ n = p„/, n /(l 

Pn' ,n' ) • 

The suppression factor S n is given by 



S n — 



dv(p„ n (a n ,v — v n )e 



(7) 



where tp Vn (a n , v — v n ) — (Avjj,n)~ 1 4'( x ) is the Voigt profile 
with x — (y — Vn) I Ai/o,n and a n is the ratio between the 
damping width and Doppler width for Lyn photons. The 
optical depth is given by 



t„{v, z;z a ) = oy 



dz' -^-n H (z') 
dz' 



(8) 



a n (T'), v 



1 — ty/c 1 + z' 



1 — V r /C 1 + Z 



2.2 The halo emissivities 



The emissivities depend on the populations n n> ; of the ex- 
cited states, where the subscripts indicate the principal 
quantum number and orbital angular momentum, respec- 
tively, of a given fine-structure state of an excited hydrogen 
atom. The equations governing the populations are, for the 
ground state m, 



dni 
~~dt 



= 2Z n n> ±A nl 



.10 - ni 



71 = 2 

and for excited states, 



dn n ;i 
dt 



niP n Si tl + E E n n',i' A n 

n'=n + l l'—l±l 

n-1 

—n n ,i E E A n l,n'l'i 

n' = l i' = i±l 



(9) 



(10) 
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where 5iji = 1 for I = I' and vanishes otherwise. Here, 
A n i n n' is the spontaneous decay rate from level n, I to level 
n'J'. 

The level populations will establish their steady-state 
values on the timescales P^ 1 ■ The system of steady-state 
equations is solved by 



n n ,i = 22 mP n i 



En— X yi 4 



(11) 



where C n n',ni is the cascade matrix expressing the probabil- 
ity that an upper state n' , I' cascades down to a lower state 
n, I via all possible transition routes. It is given by 



C„'l, n l — y ] y ] C n 1 \ ,ri 1 1" Qn" l" 



(12) 



n"=n+l i" = l±l 



with C nVM 

m' — 1 



= 5i : ii, where a n , VM = 

I X/n"=i Si"=i'±i A n >i'.n"i " is the bra nching 
ratio from level ra'Z' to level nZ (cf. lOsterbrockl (Il989h) . Th e 
transition rates A n >ii tn i are computed as in iMeiksinl (|2010h . 

The resulting emissivity due to all transitions from n to 
n' with mean frequency v nn i is 

n ^-l A nl,n'V (13) 



47T 

4-k 



I l'=l±l 



-Til 



where p n i = X^n^n Pn'C n >i,ni- The observed specific inten- 
sity is then 



Jv 



dz ~dz ^ + ^ 3 e ™ n '^ p ) 

X< AW (l + — ) (1 + Z)- ^n'] 



(14) 



for emission along the line-of-sight path l p with line-of-sight 
velocity v\\ . 



3 INFRA-RED EMISSION TOMOGRAPHY 

The optical depth is sensitive to the density, temperature 
and peculiar velocity fields. As a consequence, fluctuations 
in these quantities between the source and a given gas par- 
cel will produce fluctuations in the scattering rate through 
Eqs. @ and ([5]). These in turn will produce fluctuations in 
the emissivity of the resulting cascade radiation, opening up 
the opportunity to measure the density, temperature and 
peculiar velocity fields on small scales surrounding a source. 



3.1 Numerical simulation 

In order to estimate the magnitude of the effect of the fluc- 
tuations on the received intensity, we use a numerical cos- 
mological ACDM simulation of the IGM to compute the 
signal predicted from the surroundings of a source placed in 
the simulation volume. The simulation was run using ver- 
sion 1.0.1 of ENZC0 with the cosmological parameter values 
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Figure 1. (Clockwise from the top left) The temperature, tangen- 
tial (single component) peculiar velocity, radial peculiar velocity 
and density fields of IGM around a radiating source at the origin. 
(Distance units are proper.) 



Q, m = 0.236, il v = 0.764 and Q b = 0.041 for the matter, vac- 
uum energy and baryonic matter density parameters, a Hub- 
ble constant Ho = 100/tkm s -1 Mpc - with h = 0.73, and 
the initial BBKS power spectrum with n = 0.95 normalised 
to as = 0.74, generally consistent with CMB fluctuations 
as measured by the Wilkinson Microwave Anisotropy Probe 
|Larson et al.ll201ll ). The simulation was run in a volume 
10/i -1 Mpc on a side with 256 3 cells and 512 3 dark matter 
particles. The volume was evolved from an initial redshift 
of Zi — 530 to z = 10. The initial mean temperature is set 

to 550 K (j 1 ^) ^ 3800 K, and evolves to T ~ 2K by 
z = 10. Radiative cooling permits collapse of gas in halos, 
but there is no feedback from star formation. 

A source was assigned to the highest gas density peak 
within the simulation volume at z — 10, as recorded on the 
gas grid with cells of size 10/256 Mpc hT 1 ~ 40kpc ft -1 per 
side comoving. The corresponding density, temperature and 
peculiar velocity fields in a quadrant around the source are 
shown for a slice in Fig. \T\ with representative line-of-sight 
values shown in Fig. [2] 

The source is modelled as a starburst galaxy with a 
specific intensity approximated as fiat at the relevant fre- 
quencies, so that the source intensity a luminosity distance 
tl away is 



(4tt 



t 2~2 



(15) 



where L VL is the luminosity at the Lyman edge frequency 
vl- A continuous star formation rate lOMoyr -1 is assumed, 
with a Salpeter initial mass function for the range 1 < 
M < IOOMq and a metallicity Z — O.ObZ^ , corresponding 



to L v r = 3.8 x 10 21 Js _1 Hz -1 after 10 7 yr (|Leitherer et all 
1 1999(1 . 
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Figure 2. (Clockwise from the top left) Representative line-of- 
sight variations in the temperature, peculiar velocity, total veloc- 
ity (including Hubble flow) and density of IGM around a radiating 
source at the origin. The dotted lines are the corresponding val- 
ues for a homogeneous expanding medium. (Distance units are 
proper.) 
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Figure 4. The relative difference between the Hn emissivity along 
the ray in Fig. [2] and the corresponding rate assuming a homoge- 
neous expanding medium. 
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Figure 3. The relative difference between the Lyn scattering 
rate along the ray in Fig. [2] and the corresponding rate assuming 
a homogeneous expanding medium. 



3.2 Emission line haloes 

The resulting scattering rates of Lyn photons differ substan- 
tially from those for a homogeneous expanding medium, as 
shown in Fig. [3] The contributions due to the density, tem- 
perature and velocity departures from the mean are broken 
down in the figure. While the differences are primarily due 
to the density fluctuations, fluctuations in the temperature 
and expansion velocity contribute non-negligibly. In particu- 
lar, the scattering rate becomes increasingly sensitive to the 
peculiar velocity field towards the higher orders, for which 
scattering in the Doppler core dominates increasingly over 
Lorentz wing scattering. As the optical depth increases with 
decreasing total velocity gradient, regions with a peculiar 
velocity gradient dv/dr < tend to have a suppressed scat- 
tering rate compared with regions with dv/dr > 0. This may 
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Figure 5. Hn emissivity along the ray in Fig. [2] 



be seen by comparing the fluctuations in the scattering rate 
in Fig. [3] with the slope in the total velocity shown in the 
lower right panel of Fig. [2] 

The Balmer emissivity profiles for Ha through H<5 are 
shown in Fig. |4]for the ray shown in Fig. [2] The fluctuations 
follow those in the Lyn scattering rates, but are somewhat 
suppressed in magnitude, as shown in Fig. [5] This is be- 
cause, while the emissivity is proportional to the local gas 
density, so that fluctuations in the emissivity include a linear 
dependence on the density fluctuations, the scattering rate 
fluctuates oppositely to the density fluctuations because of 
the effect of the optical depth variations. The sum of the 
contributions to the relative fluctuations in the emissivity is 
thus reduced compared with the magnitude of the relative 
density fluctuations. 

The derivation of the surface brightness is computa- 
tionally intensive, as the cascade equations must be solved 
separately within each simulation cell. For this reason the 
computation is restricted to a plane perpendicular to the 
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Figure 6. Differential Balmer surface brightness fluctuations 
compared with the homogeneous expanding IGM case allowing 
only for the density fluctuation contribution. 
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Figure 7. Differential Balmer surface brightness fluctuations 
compared with the homogeneous expanding IGM case allowing 
only for the density and temperature fluctuation contributions. 



line-of-sight. A narrow-band filter would include contribu- 
tions from the nearby foreground and background regions 
as well. Modelling a specific filter arrangement would re- 
quire further layers to be included for a realistic estimate, 
which would tend to blur the image if the filter width corre- 
sponded to a length scale that exceeds the coherence scale of 
the inhomogeneities. The results shown here are thus only 
representative of the magnitude of the effects that should 
arise from inhomogeneities in the density, temperature and 
peculiar velocity fields. 
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Figure 8. Differential Balmer surface brightness fluctuations 
compared with the homogeneous expanding IGM case allowing 
only for the density and velocity fluctuation contributions. 
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Figure 9. Differential Balmer surface brightness fluctuations 
compared with the homogeneous expanding IGM case allowing 
for the effects of density, temperature and velocity fluctuations. 



The effects of the density, temperature and velocity fluc- 
tuations on the Balmer surface brightnesses are illustrated 
in the differential maps shown in Figs. [6] [7] [8] and [9] The 
maps show the differences from the case for a homogeneous 
expanding IGM. The differential map allowing only for the 
density fluctuations shows the suppression by the increased 
optical depth in dense regions, as in Fig. [5] Adding in the 
contribution of the temperature fluctuations substantially 
enhances the structures. Adding the peculiar velocity con- 
tribution to the density emphasizes differentials in the struc- 
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Figure 10. Balmer surface brightnesses including the effects of 
density, temperature and velocity fluctuations. 



Source 


Hcv 
5.91^tm 


H/3 
4.38/im 


H7 
3.91/im 


H<5 
3.69/xm 


Zodiacal 


4.7 x 10 4 


7.7 x 10 3 


3.5 x 10 3 


2.6 x 10 3 


Halo 


0.048 


0.018 


0.011 


0.0089 


SNR/Mioo 


0.16 *4 /2 


0.15 t\ /2 


0.14 t 1 / 2 


0.13 t\ /2 



Table 1. Estimated photon flux (photons m~ 2 s —1 ) from low- 
level zodiacal light and the integrated halo flux within one ar- 
cmin of a starburst at Z = 8, through an aperture of solid an- 
gle 4 arcmin 2 and narrow- band filter with resolution A/AA = 
100. The starburst is normalised to a star formation rate of 
IOOM100 M yr" 1 . The last row provides the signal-to-noise ratio 
after an integration time (4 in units of 10 4 s for flux incident on a 
mirror and detector with effective collecting area 50 m 2 . An IGM 
temperature T = 10 K is assumed. 



tures, allowing the possibility of tracing the peculiar velocity 
field around forming galaxies on small scales. 

The combined map is shown in Fig. 1101 The rippling 
effect of the IGM is clearly discernable. The influence of the 
peculiar velocity becomes increasingly strong for the higher 
orders, where the dominant scattering producing the signals 
moves from the Voigt wings to the Doppler core. 



4 DETECTION OF INFRA-RED EMISSION 
HALOES 

The Balmer emission lines produced in the diffuse IGM by 
galaxies during the Epoch of Reionization will redshift to 
the infrared. In this section, an estimate is made for the 
requirements to detect the haloes. It is assumed the detector 
is limited by photon noise, which will be dominated by the 
IR background, primarily zodiacal light for a telescope in 
space. A fiducial collecting area of And = 50 m 2 is adopted 



45 £ 

-16 <3 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 
x [Mpc'l 



with a field of view of FOVfld = 2x2 arcmin 2 . Imaging 
using a narrow-band filter of resolution R — A/AA = 100 is 
assumed. 

Estimat es for the zodi acal light background^ based on 
the model of IWrightl (|l998h are converted to count rates in 
Table [T] The intensity of zodiacal light depends on the posi- 
tion and pointing direction of a satellite. The lower estimates 
17 g are quoted here. 

Lyman Break Galaxies and Lya emitters at z > 
6 suggest star formation rates of up to 30 M yr _1 
-n l|Gonzalez et all l201ll ; iMcLure et"afl l201lh . Much higher 
i3_„ star formation rates are known in the lower redshift Uni- 
|4 r, s verse, reaching values of 500 — 1000 Mq yr _1 , but these are 

15 g: associated with ultra-luminous infra-red galaxies (ULIRGs) 

16 0? in dusty environ ments, which sev erely attenuate the ultra- 

17 § violet radiation (|Kennicuttl 1 199&t ) . On the other hand, by 
[ z — 1 — 2 ULIRGs are found to be much more transparent 

to UV radiation (Dad dTet al.ll2007h . Accordingly, our esti- 
mates are normalized to a source with a star formation rate 
of 100 M Q yr" 1 . 

The resulting lower order Balme r fluxes a re sho wn in 
Table [__ based on the estimates in iMeiksinl l|2010l ). The 
corresponding signal-to-noise ratios (SNR) after integrating 
over a time t are provided in the last row of the table. The 
SNR is nearly independent of the Balmer order. Increasing 
the resolution to R = 500 would be adequate for capturing 
the emitting flux while minimizing the zodiacal light back- 
ground. A starburst as great as 700 Mq yr _1 would still be 
required to achieve a 3er detection in 2 x 10 4 s. Allowing for 
the stacking of multiple bands in an instrument that split 
the beam into a range of bands tuned to correspond to the 
different Balmer series, would reduce the required star for- 
mation rate for detection by the square-root of the number 
of bands stacked. 

The haloes could plausibly be detected by the James 
Webb Space Telescope (JWST) using the Tunable Filter Im- 
ager for a sufficiently bright source. The narrowest filter 
width, with resolutioijf] R ~ 100, is sub-optimal, as it is 
broader than the emission feature, and thus lets in an exces- 
sive amount of background light compared with the signal. 
The wavelength range of the detector would also miss Ha 
for a source at redshift zs > 6.6, however it would capture 
the higher orders. The mirror has a collecting area of 25 m 2 . 
Allowing for a photon throughput of 50 per cent., a 5a de- 
tection of H/3 for a 1000 Mq yr" 1 source would require an 
integration time of 4.7 x 10 s. 



5 CONCLUSIONS 

The search for the highest redshift galaxies may have identi- 
fied forming galaxies responsible at least for a large compo- 
nent of the reionization of the intergalactic medium. Estab- 
lishing that the IGM was actually neutral, however, requires 
direct detection of the neutral medium. A major goal of ex- 
isting and planned radio facilities is to detect the neutral 
IGM 21cm signal from the Epoch of Reionization. A com- 
plementary path is to discover the UV light from the sources 



At http://jwstetc.stsci.edu/ctcstatic/users_guide/t8_background.html 

At http:// www.stsci.edu/jwst/doc-archive/handbooks/JWST_Primer_v20.pdf 
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re-radiated by the IGM and received as infra-red radiation. 
Rescattered Lya photons offer one possibility, although the 
emission is restricted by the spatial diffusion of the photons 
to the vicinity of the source. Higher energy photons will 
redshift into local higher order Lyman series photons, excit- 
ing extended Balmer and higher order emission line haloes 
through electron cascades following Lyman photon scatter- 
ing. 

We have shown that the fluctuations in the signal are 
sensitive not only to the density fluctuations in the sur- 
rounding gas through the mean free path of the scattered 
Lyman photons, but of the velocity field as well. This offers 
the possibility of mapping out both the ionization struc- 
ture and the peculiar velocity field, produced by inflows or 
wind-driven outflows, surrounding the earliest major radia- 
tion sources, whether galaxies or AGN. 

The principle obstacle to the detection of the infra-red 
haloes is the zodiacal light background. The ideal observing 
instrument would be a space telescope with an effective col- 
lecting area of 50 m 2 , a field of view of a few to several 
square arcminutes, and a detector sensitive to the wave- 
length range 1-7 pm with the capability of simultaneous 
imaging in several narrow tunable bands with a resolution 
of R ~ 500. Basing an estimate on a pointing to a region of 
sky with a low-level of zodiacal light background, we find a 
700 Moyr -1 starburst galaxy could be detected at the 3a 
level in a single band in a 2 x 10 4 s integration, or a star- 
burst of half the strength if observed simultaneously in four 
bands. 

Although sub-optimal in design for this purpose, the 
Tunable Filter Imager on JWST could detect a z — 8 
1000 M Q yr _1 starburst galaxy at the 5a level in H/3 in a 
5 x 10 5 s integration. While a major resource investment, 
it is half the time allocated to the Hubble Space Telescope 
Ultra-Deep Field and would offer the prize of a definitive 
detection of the Epoch of Reionization. 
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